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We present a detailed infrared study of the insulator-to-metal transition (IMT) in vanadium diox- 
ide (VO2) thin films. Conventional infrared spectroscopy was employed to investigate the IMT in 
the far-field. Scanning near-field infrared microscopy directly revealed the percolative IMT with in- 
creasing temperature. We confirmed that the phase transition is also percolative with cooling across 
the IMT. We present extensive near-field infrared images of phase coexistence in the IMT regime in 
VO2. We find that the coexisting insulating and metallic regions at a fixed temperature are static 
on the time scale of our measurements. A novel approach for analyzing the far-field and near-field 
infrared data within the Bruggeman effective medium theory was employed to extract the optical 
constants of the incipient metallic puddles at the onset of the IMT. We found divergent effective 
carrier mass in the metallic puddles that demonstrates the importance of electronic correlations to 
the IMT in VO2. We employ the extended dipole model for a quantitative analysis of the observed 
near-field infrared amplitude contrast and compare the results with those obtained with the basic 
dipole model. 



PACS numbers: 71.30.+h, 71.27.+a, 78.20.-e, 78.30.-j 

I. INTRODUCTION 

Vanadium dioxide (VO2) undergoes an insulator-to- 
metal transition (IMT) at T c « 340 K. The IMT is accom- 
panied by a change of lattice structure from a monoclinic 
(Ml) unit cell in the low-temperature insulating phase to 
tetragonal (rutile) symmetry in the metallic phase. VO2 
has been studied extensively for two main reasons. First, 
the driving mechanism of the IMT has remained con- 
troversial and many experimental and theoretical studies 
have focused on understanding it Pj^j S econdly, VO2 has 
immense potential for applications^ 4 ^ because the phase 
transition occurs only 40 K above room temperature and 
involves dramatic changes in transport, infrared and op- 
tical properties. Moreover, the phase transition can also 
be triggered optically^' and electrically^^ leading to new 
optical and electronic devices. 

The reason the driving mechanism of the IMT in 
VO2 has been a subject of controversy is the complex 
nature of the interactions in this material. VO2 con- 
sists of just two elements, vanadium and oxygen, with 
well-documented lattice structures in the Ml and rutile 
phases. The model of delocalized rf-electrons of vana- 
dium (one per vanadium ion) moving in the periodic po- 
tential created by the lattice was expected to explain the 
documented features of the IMT in VO2. In the early 
literature, the phase transition in VO2 was described in 
terms of the structural distortion that leads to the for- 



mation of vanadium pairs and unit-cell doublingP This 
Peierls instability was thought to occur due to the cou- 
pling of th e condu ction electrons to a soft phonon mode of 
vanadium! 10 * 11 * 12 ^ However, the mod el of a Peierls insta- 
bility within band theory alone cannot satisfactorily 
explain several experimental aspects of the phase tran- 
sition and the properties of the insulating and metallic 
phases existing in its vicinity. Following are a few ex- 
amples of such experimental observations: the » 0.6 eV 
energy gap in the Ml insul ator measur ed by infrared 
and photoemission technique d 15 * 16 * 17 * 18 * 19 * ; the changes in 
optical conductivity a nd spectral weight up to at least 
6 eV due to the iM fl 15 ' 16 ' 17 ' 1 ^ the satellite feature 
at w 1 eV below the Fermi energy seen in photoemis- 
sion in the rutile metal which may well be a remnant of 
the lower Hubbard bancP^ and lack of resistivity satura- 
tion in the rutile metal beyond the Ioffe-Regel-Mott limit 
of metallic transport which is difficult to explain within 
a framework based on electron-phonon coupling ! 18 * 21 * 22 * 
There is growing evidence that electron-electron corre- 
lations cannot be neglected in any viable description of 
the IMT in vOz! 18 ' 19 ' 23 ' 24 *^ 1 ™ Several previous the- 
oretical stu dies have em phasized the role of electronic 
correlations.^ 8 ' 29 ' 30 ' 31 ' 32 ' 

It appears that the interplay among various electronic 
degrees of freedom (charge, orbital and spin) and their 
coupling to the lattice degrees of freedom (phonons) 
leads to complex behavior in V02- This is a hallmark 
of many transition-metal oxides that exhibit remark- 
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able properties like high-temperature superconductivity, 
colossal magnetoresistance and multiferroicityP^ Many 
of these transition metal oxides have a rich phase dia- 
gram with a diversity of phasesPS Some of the exotic 
properties may well be due to phase coexistence close to 
phase boundaries. Recent experimental work on high- 
temperature superconductors and manganites provides 
direct evidence of phase separation and inhomogcneity 
in the se materials on mesoscopic and nanoscopic length 
scales. 13 ^ 1 - 36 1 37 1 38 1 39 1 It is thought that competing elec- 
tronic and phononic degrees of freedom lead to phase sep- 
aration and a multiplicity of phases in transition metal 
oxides The evolution of one phase to another appears 
to proceed on the nanometer scale leading to phase coex- 
istence and inhomogeneity. This motivated us to inves- 
tigate the IMT regime in VO2 on the nanometer length 
scales with the pioneering technique of scattering scan- 
ning near-field infrared microscopy (s-SNIM) 

While the IMT in VO2 has been stu died previously 
by far-field infrared and THz technique d 41 * 42 !, our stud- 
ies combine far-field infrared spectroscopy with near-field 
infrared imaging!^ With s-SNIM, we have imaged the 
evolution of the thermally-induced IMT in a thin film 
of VO2. We found that the VO2 film underwent a per- 
colative transition with increasing temperature across the 
IMT. We analyzed the s-SNIM data along with the far- 
field infrared spectroscopy data within the Bruggeman 
effective medium theory (EMT) , and discovered that the 
metallic islands embedded in the insulating host exhibit 
divergent effective mass which is evidence of th e Mott 
transition within the Brinkman-Rice framework! 23 ! 43 * 4 ^ 
In this paper, we further discuss the underlying physics 
of the IMT in V0 2 . We present extensive s-SNIM im- 
ages obtained while heating and cooling across the IMT 
which confirm the reproducibility of our results. We 
demonstrate the static nature of the coexisting insulat- 
ing and metallic regions over many minutes at a fixed 
temperature. A quantitative discussion of the observed 
near-field infrared amplitude contrast in the IMT regime 
is included. The amplitude contrast was modeled with 
the so-called dipole description^ and its refined version. 
Moreover, we explore the novel approach to data analy- 
sis of the far-field and near-field infrared data within the 
Bruggeman effective medium theory and expound upon 
the merits and limitations of this procedure. 

This article is divided into several sections. The details 
of the sample preparation and characterization are given 
in the next section. The far-field infrared spectroscopy 
experiments and data analysis are described in Section 
III. This is followed by a description of the s-SNIM tech- 
nique and presentation of the data in Sections IV and 
V respectively. The data analysis within the Bruggeman 
effective medium theory is described in Section VI. This 
is followed by a Summary and Outlook Section and three 
appendices. 



II. SAMPLES 

VO2 films about 100 nm thick were grown on (1012) 
oriented sapphire (AI2O3) substrates by the sol gel 
method. The details of growth and characterization are 
given in Ref. 135] In Fig. [TJi we present the X-ray diffrac- 
tion data on the VO2 film which shows that the film is 
single phase MI-VO2 with (200) orientation. The resis- 
tance of the film was measured by a standard four-probe 
method, and the data is displayed in Fig.[T|3. The film un- 
dergoes a first-order, hysteretic IMT with four orders of 
magnitude change in the resistance. We chose thin films 
of VO2 for infrared measurements in the phase transi- 
tion regime instead of single crystals because the IMT is 
destructive for VO2 crystals leading to potential prob- 
lems of reproducibility and significant uncertainties in the 
quantitative results in the phase transition regime. Data 
obtained for high quality films grown on lattice matched 
substrates are free of the above complications and thus 
enable studies of the intrinsic properties of V02- Our 
VO2 films do not show signs of deterioration even af- 
ter going through several cycles across the IMT. This 
was verified by repeated resistance measurements of the 
IMT, and by ellipsometric and reflectance experiments 
that give reproducible data in the insulating and metal- 
lic phases before and after cycling through the IMT. 



III. FAR-FIELD INFRARED SPECTROSCOPY 

Spectroscopic ellipsometry provides the ellipsomet- 
ric coefficients {^(w) and A(cj)) at each measured 
frequency.^ This enables us to obtain the optical con- 
stants of the VO2 film without recourse to Kramers- 
Kronig analysis which may not be applicable in the IMT 
regime of VO2 because of inhomogeneity and phase- 
coexistence. Provided the incident wavelength is large 
compared to the size of the inhomogeneities, the inho- 
mogeneous system can be described by an effective di- 
electric function 4 ^. This is certainly expected to be the 
case for long wavelengths in the far- and mid-infrared 
regime (between 40 cm -1 and 5000 cm^ 1 ). Further- 
more, we find that in accord with the EMT^I, weighted 
averaging of the dielectric functions of the insulating 
phase and rutile metallic phase at shorter wavelengths 
(higher frequencies) adequately describes the data (see 
Appendix A). Therefore, it is appropriate to assign an 
effective dielectric function to the inhomogeneous, phase- 
separated regime of VO2 even when the wavelength be- 
comes comparable to the size of the metallic puddles. 
The effective dielectric function of the inhomogeneous, 
phase-separated regime of VO2 is thus obtained via the 
standard analysis of the combined ellipsometric and re- 
flectance data based on a two-l ayer model of a VO2 
film on a sapphire substrate! 19 * 47 ! The ellipsometric data 
spanned the frequency range 400 - 20000 cm -1 , and near- 
normal incidence reflectance covered the frequency range 
40 - 680 cm -1 . 
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We note that in general, the complex dielectric func- 
tion i(uS)= ei(tL>) + j£2(w) of a material is related to the 
complex conductivity d(ui)= <J\(uj) + icr 2 (o;) by the fol- 
lowing equation: 



i(w) = 1 + er(«) (1) 

The real part of the optical conductivity <j\ (uj) is plot- 
ted in Fig. as a function of frequency (oS) for repre- 
sentative temperatures in the IMT regime of VC>2- The 
inset in Fig. [TJ3 shows the temperature-dependence of the 
real part of the dielectric function e\(u> = 50 cm -1 ) in 
the IMT regime. Fig. [TJi and its inset show <j\(oj) and 
ei(o; = 50 cm^ 1 ) respectively with phonon contributions 
subtracted. The contributions of VO2 phonons to the di- 
electric functions were modeled by Lorentzian oscillators 
and subsequently subtracted so that the electronic contri- 
bution to the optical constants could be unambiguously 
presented and analyzed. 

An optical gap of w 4000 cm -1 (0.5 cV) is evident 
in <ti(u>) of the insulating phase (295 K < T < 341 
K). As the V0 2 film enters the IMT regime for T > 
341 K, the conductivity in the gap increases and grad- 
ually fills the gap. The first inter-band transition cen- 
tered at 10500 cm -1 does not appear to shift to lower 
frequencies as would be expected if the gap were to col- 
lapse. The filling up of the optical gap, rather than a 
collapse, is also seen in other correlated (or Mott) in- 
sulators like the cuprate s that are rendered conducting 
through chemical doping P 49 * 50 * In addition, we observe 
an isosbestic point in cr 1 (w) characterized by invariant 
conductivity (depicted by an open circle in Fig. [lj:,d) 
The isosbestic point is also seen in other Mott systems^ 
although its significance is yet to be understood. The 
observations discussed in the preceding sentences point 
to the importance of electronic correlations in VO2. As 
the temperature is increased, the low frequency conduc- 
tivity increases until the film becomes fully metallic as 
seen from the broad Drude-like feature at T = 360 K. 
We note that the enhanced spectral weight at low ener- 
gies as the film becomes increasingly metallic is borrowed 
from higher lying inter-band transitions! 17 * 18 * 1 ^ The di- 
vergence of ei shown in the insets of Fig. |Tj3,cl is similar 
to that known from percolative insulator-to-metal transi- 
tions where it occurs due to enhanced capacitive coupling 
between the metallic region s as they grow and proliferate 
in the insulating hostP^^ ei abruptly becomes negative 
at T = 343 K signaling that the metallic regions have 
percolated (connected) forming a macroscopic conduct- 
ing path for the charge carriers. 

IV. SCANNING NEAR-FIELD INFRARED 
MICROSCOPY 

The characteristic evolution of the effective area- 
averaged optical constants in the IMT regime of V0 2 



provides indirect evidence of a percolative phase transi- 
tion with coexisting insulating and metallic regions. Our 
objective was to directly observe the percolation with 
scanning near-field optical techniques. We chose to di- 
rectly image the IMT with s-SNIM at the mid-infrared 
frequency of 930 cm -1 . As can be seen from Fig. [lj;,d, 
the optical conductivity of the insulating and metallic 
phases of V0 2 is very different at this frequency and so 
a large contrast is expected in the near-field amplitude 
signal. 

Scattering scanning near-field infrared microscopy was 
performed with a custom Atomic Force Microscope 
(AFM) using Pt-coated silicon tips in the tapping 
modeP2lSni The u> = 930 cm -1 infrared frequency was 
provided by a CO2 laser. The pseudo-heterodyne scheme 
with second harmonic demodulation was employed to iso- 
late and detect the near-field amplitude signal.^ The 
spatial resolution of the near-field infrared probe is typ- 
ically 10-20 nm and is set by the radius of curvature of 
the tip. The near-field penetrates about 20 nm into the 
VO2 film and can therefore be considered a bulk probe 
of the insulator-to-metal transition!^] 

The near-field interaction between the tip and the sam- 
ple is described by the dipole model. This is the simplest 
model that quantitatively describes the near-field inter- 
action and has been found to be accurate in some cases 
Within this model, the complex scattered signal at the 
second harmonic of the tapping frequency £1 is a function 
of the optical constants of the tip and the sample.^ The 
incident infrared light polarizes the tip and within the 
dipole model, it is sufficient to model the polarizability 
of the tip with an effective spherical dipole of radius a 
(« 20 nm). The interaction between the tip dipole and 
the sample is modeled by considering a mirror dipole in 
the sample whose polarizability depends on the local op- 
tical constants of the sample in the region directly below 
the tip. The effective polarizability of the interacting 
tip-sample system (a e //) is given by: 

1 167r(a+z(t))3 

a e f f is complex- valued and relates the incident Ei and 
scattered E s electric fields via E s cx a e ffEi. The polar- 
izability of the tip a is given by a — 4ira 3 (i t — l)/(e t + 2) 
where it is the complex dielectric function of the tip 
at the incident laser frequency. The polarizability of 
the mirror dipole in the sample is given by a(3 where 
(3 = (i s — 1) /(i s + 1) expresses the local dielectric re- 
sponse of the sample that depends on its complex dielec- 
tric function e s at the incident laser frequency. 

The effective polarizability of the tip-sample system 
is a non-linear function of the tip-sample distance z. 
The tip-sample distance varies sinusoidally with time 
z(t) = Zo(l + cos(f2i)) due to oscillation of the tip. In the 
tapping mode, the tip oscillates with an amplitude z of 
about 40 nm at a frequency SI = 25 kHz. Therefore, the 
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effective polarizability and hence the near-field signal is 
modulated at harmonics of f2. 



a e f f = (a e //) + (a eff ) 1 cos(flt) + (a e //) 2 cos(2fii) 
+(a e //) 3 cos(3f2i) + ... + (a e //) ri cos(nf2t). 

(3) 

The nth harmonic coefficient of the near-field polariz- 
ability is as follows: 



(a 



eft) 



n = — / a e // cos (nfit)di. 

27r J-tt/O 



(4) 



Here a e // is given by eq. [2] The nth harmonic of 
the scattered electric field (Eg) n = s„e ! *" oc (a e ff) n Ei, 
where s n and 0„ are the nth harmonic amplitude and 
phase coefficients respectively. Thus the second harmonic 
coefficient of the near-field amplitude S2 can be calculated 
(up to a proportionality constant that can be calibrated) 
provided the complex dielectric functions of the tip and 
sample at the incident laser frequency are known. 

More recently, it has been suggested that the (1+/3) 
term in the numerator of eq. [2] be replaced by the term 



(1 



coefficient^^ 



where r is the far-field Fresnel reflectance 



a eff 



a/3 

167r(a+z(t)) 3 



(5) 



In the above equation, the (1 + r) 2 term accounts for 
the indirect illumination of the tip from reflections at the 
sample surface. 

The dielectric function it of Pt at ui = 930 cm -1 is 
-1300 + 960iP The dielectric functions of V0 2 at w 
— 930 cm' 1 are obtained from ellipsometric chracteri- 
zation. These are as follows: 4.9 with negligibly small 
imaginary part for the monoclinic insulating phase, and 
—35 + 119z for the rutile metallic phase. Therefore, the 
scattered near-field amplitude in the metallic phase of 
VO2 is expected to be about twice that in the insulating 
phase from basic dipole model (eq. [2]). However, the ex- 
tended dipole model (eq. [5| gives a factor of 3 ratio for 
the near-field amplitudes in the metallic and insulating 
phases. 



V. INFRARED NANO-IM AGING OF THE 
INSULATOR-TO-METAL TRANSITION IN VO2 

We now turn to Fig. [2] which displays the second har- 
monic amplitude images in the 1MT regime of VO2 ob- 
tained with s-SNIM at the incident laser frequency of 
930 cm -1 . At T = 34f K, we observe an essentially 
uniform amplitude signal. As the temperature is in- 
creased, nanometer-scale regions with higher near-field 
amplitude appear. The increased near-field amplitude 



of these nano-scale puddles is due to their metallicity. 
The metallic puddles subsequently proliferate and grow 
in size, and then coalesce with increasing temperature. 
The T — 360 K near-field map shows higher near-field 
amplitude indicating that the scanned region has be- 
come completely metallic. Near-field amplitude images 
were also obtained while cooling the VO2 film through 
the 1MT and are shown in Appendix B. We note here 
that coexisting insulating and metallic domains in VO2 
nano-rods in the IMT regime have been observed with 
optical microscopy in Ref. [BOl Moreover, the IMT at the 
surface of VO2 films has also been imaged with scanning 
tunneling microscopy.^ 

The distribution of the scattered near-field amplitude 
for various temperatures in the IMT regime of VO2 is 
plotted in Fig. 3 1. This distribution is obtained from the 
images in Fig. 2 The asymmetry of the histograms and 
the presence of two peak structures are evidence of coex- 
isting insulating and metallic phases in the IMT regime 
of V0 2 . The peak at lower near- field amplitudes cor- 
responds to the insulating phase and the one at higher 
amplitudes corresponds to the metallic phase. The me- 
dian amplitude for the metallic phase at T = 360 K is 
about 15 times higher than the median amplitude for the 
insulating phase. This is far higher than the factor of 2-3 
expected from the dipole model. While eq. [5] improves 
the description of the near-field contrast in our data com- 
pared to eq. [2j it still underestimates the contrast com- 
pared to experiment. A similar discrepancy between the 
measured near-field contrast and that predicted by the 
dipole model has also been noted for near-field terahertz 
measurements on the IMT in V0 2 P^1 

We note that each of separate images in Fig. [2] has 
only a factor of 2-3 variation of the near-field amplitude 
contrast. Moreover, the median position of the "insulat- 
ing" peak in Fig. [3^, shows a small but noticeable increase 
with temperature without apparent reason. These obser- 
vations led us to normalize the measured amplitude to 
the peak position of the insulating phase. The resulting 
distributions plotted as a function of the normalized am- 
plitude are in fairly good agreement with the predictions 
of the dipole models (see Fig. [3}d). The amplitude from 
the metallic phase is now a factor of 2-4 higher than that 
from the insulating phase. Note that the basic dipole 
model (eq. [2]) gives better quantitative agreement in the 
regime with smaller sizes of the metallic puddles in the in- 
sulating host and vice versa. The extended dipole model 
(eq. [5]) explains the data well in the regime where the in- 
sulating and metallic regions have similar length scales. 
The apparent continuous increase in the near-field ampli- 
tude of the metallic phase with increasing temperature 
(see Fig. [3^,) and the 15 times higher near-field ampli- 
tude from the metallic phase compared to the insulating 
phase are phenomena that are not fully understood at 
present and require further investigation. The insulat- 
ing and metallic fractions can, nevertheless be extracted 
from the near-field images and amplitude distributions 
as explained in the next section. 
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We note that the near-field images are repeatable when 
taken in five minute sequences over the same sample area 
and at a fixed temperature because they show nearly 
identical patterns of metallic puddles in the insulating 
host (see Appendix B). This indicates that the possible 
effect of dynamic fluctuations on the static patterns is 
small. This is discussed in detail in Appendix B. 



VI. BROAD-BAND INFRARED RESPONSE OF 
THE METALLIC REGIONS IN THE PHASE 
COEXISTENCE REGIME OF V0 2 



T (K) / 



342 0.18 0.2 
342.6 0.31 0.33 

343 0.48 0.45 
343.6 0.7 0.5 



TABLE I: The values of the filling fraction / and depolar- 
ization factor q used to obtain the optical constants of the 
metallic regions at selected temperatures within the frame- 
work of the EMT formalism. 



As shown in Appendix A, weighted averaging of op- 
tical constants of the insulating phase and of the rutile 
metallic phase within the Bruggeman effective medium 
theory (EMT) produces a good description of the far- 
field infrared data near the onset of the insulator-to- metal 
transition at high frequencies. Deviations from this de- 
scription become progressively stronger at longer wave- 
lengths, a regime where the EMT formalism is expected 
to be best applicable. 48 This implies that the optical con- 
stants of the metallic puddles when they first appear at 
the onset of the insulator-to-metal transition are differ- 
ent from the optical constants of the rutile metal at T = 
360 K. Therefore, we use the EMT to extract the optical 
constants of the metallic puddles. 

Within the Bruggeman effective medium theory,^ the 
effective (complex) optical constants e E {oj) of a two- 
component inhomogencous system are given by: 



i a (w) - e~ E (uj) . . e b (u>) - e~ E (uj) 



£ b (cj) + ^e E (u;) 



In this equation, e a (u>) and £b(oj) are the complex opti- 
cal constants of the metallic and insulating phases respec- 
tively, / and (1-/) are the volume fractions of the metallic 
and insulating phases respectively, and q is the depo- 
larization factor that depends on the shape of the com- 
ponents which is inferred from near-field images. This 
factor is taken to be 0.2-0.4 assuming nearly spherical 
metallic regions at low concentrations of the newborn 
metallic state. This is because their lateral dimensions 
are less than the film thickness and the out-of-plane di- 
mension is assumed to be similar to the lateral dimen- 
sions. The depolarization factor continuously increases 
to 0.5 (appropriate for thin flat disks) at higher concen- 
trations of metallic clusters as their lateral dimensions 
exceed the out-of-plane dimension which is limited by 
the thickness of the film. The volume fractions of the 
insulating and metallic phases are obtained directly from 
the near-field images by fitting the amplitude distribu- 
tion curves in Fig. [3^ to two Gaussians, with the area 
under each Gaussian yielding the volume fraction of the 
respective phase. The optical constants of the insulat- 
ing phase e b (u)) and the inhomogeneous transition regime 
e E {uj) have also been measured (see Fig. 1). The optical 



constants of the insulating phase are assumed to be in- 
dependent of temperature in this analysis. Since e~ E (cj), 
and / are known whereas the range of q is con- 
strained based on the knowledge of the metallic cluster 
size and film thickness, the equation can be solved to ob- 
tain the optical constants of the metallic regions e a (uj), 
and hence the complex conductivity <7 a (uj) for each tem- 
perature. The real part of the conductivity ai a (u>) ex- 
tracted by this procedure is plotted in Fig. |2Ja) for rep- 
resentative temperatures in the IMT regime. The values 
of / and q used for the analysis are listed in Table [TJ 

We note that the form of ai a (uS) for the metallic re- 
gions that are separated by the insulating host exhibits 
a narrow Drude-like peak at low frequencies followed by 
a "dip- hump" structure at higher frequencies. This is 
quite different from ai a (uj) of the rutile metal at T = 
360 K which has a broad and featureless Drude peak. 
As the temperature is increased, and the metallic regions 
proliferate and grow, their a\ a {oj) rapidly evolves to that 
seen in the rutile metallic phase. Here we note that the 
extracted optical constants of the metallic puddles at uj 
= 930 cm -1 are different from the optical constants of 
the T — 360 K rutile metal at this frequency. Neverthe- 
less, within the dipole model (eqs|2]and[5]), the near-field 
scattering amplitude at uj = 930 cm -1 for the metallic 
puddles is expected to be nearly the same as that from 
the T — 360 K rutile metal. 

Uncertainties in the magnitude of q do affect the be- 
havior of the extracted spectra of the metallic puddles at 
the lowest frequencies accessible in our experiment. This 
is shown in Fig. [5] Changing q for a fixed value of / re- 
sults in a variation in the <J\{uj) spectra below 800 cm -1 
whereas the variation is negligibly small above 800 cm -1 . 
We see that the mid-infrared band ( "hump" ) peaked at ~ 
1800 cm -1 and the "dip" feature below this frequency are 
robust features independent of the value of q. Below 300 
cm -1 however, the choice of q produces either a Drude- 
like feature or a finite energy mode which is consistent 
with Drude dynamics modified by localization P^l The lat- 
ter behavior is often seen in systems close to the me tal- 
insulator transition in the presence of disorder. 63 6 ^We 
note that for a value of / obtained from the near-field im- 
age at a fixed temperature, the range of q values is highly 
constrained. Within EMT, percolation is predicted for / 
> g.Bal However, the far-field effective spectra and the 
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near- field patterns for T — 342.6 K indicate that per- 
colation has not yet occurred and this means that / < 
q. Since solutions of eq(6] yielding a finite energy mode 
occur for / > q, we infer that these are unlikely to be 
correct, although they cannot be dismissed entirely. We 
also see that for q > 0.34, <j\ a turns negative over a finite 
frequency range and is unphysical. Hence, the spectra 
showing a narrow Drude-like mode centered at zero fre- 
quency are the most probable candidates for the optical 
constants of the metallic puddles. 

The scattering rate 1/t(lo) and the mass enhancement 
factor m*(oj)/mi > of the metallic puddles obtained from 
the extended Drude analysis of the optical constants (see 
Appendix C) are plotted in Fig. [4jb),(c) for represen- 
tative temperatures in the IMT regime. The scattering 
rate exhibits a gap-like form at low frequencies with an 
overshoot at ~ 1000 cm -1 . We attribute the suppression 
of the scattering rate at low frequencies to partial gap- 
ping of the Fermi surface (i.e. a pseudogap). The con- 
nection between the suppression in the scattering rate 
and the opening of a pseudogap has been made earlier 
for the u nder-dop ed cuprates that are close to a Mott 
transition! 66 * 67 1 68 1 The pseudogap in VO2 could be due to 
Mott physics, the Peierls instability^ or a complex inter- 
play between the two. We cannot rule out the possibility 
that this peculiar form of \/t{lS) and the "dip- hump" 
structure in a\{uS) are due to a low-lying interband tran- 
sition. 

Note that the mass enhancement factor of the metallic 
puddles increases substantially at low frequencies. The 
mass enhancement factor in the limit of uo — > diverges 
as one approaches the IMT with varying temperature 
(see inset of Fig. Qc)). The mass divergence is likely due 
to electron-elect ron i nteractions and is a fingerprint of a 
Mott transition ! 43 * 44 ! Note that our observation of mass 
divergence is not inconsistent with the first-order nature 
of the phase transition in VO2 because about two-thirds 
of the entropy associated with the first-order transition is 
related to the lattice and the rest is electronic.^ The rela- 
tionship of electronic mass divergence and the electronic 
IMT to the structural transition needs to be explored in 
the future with nanoscale X-ray diffraction measurements 
of the metallic puddles. 

As mentioned in the preceding paragraph, the form of 
the conductivity below 300 cm -1 is somewhat ambigu- 
ous. While the narrow Drude-like solutions show a varia- 
tion in the low frequency conductivity, all such solutions 
lead to enhanced effective mass and an optical pseudo- 
gap in the metallic puddles within the extended Drude 
analysis. For solutions exhibiting a low energy mode, 
the reduced oscillator strength of this low-energy mode 
compared to that of the Drude mode in the rutile metal 
indicates enhanced optical effective mass. 

Interfacial scattering and localization effects will be- 
come important when the size of the metallic puddles is 
comparable to or less than the intrinsic mean free path 
of the charge carriers within these puddles. Our estimate 
of the mean free path based on the scattering rate data 



in Fig. |4|b) is of the order of a few nanometers. The 
smallest puddles we can detect (limited by the 20 nm 
spatial resolution of our probe) are still in the regime 
where the mean free path is shorter than the spatial ex- 
tent of metallic regions. Therefore, interfacial scattering 
and localization are unlikely to affect charge dynamics in 
the metallic puddles. However, we cannot rule out the 
possibility of quasiparticle localization which normally 
leads to a finite energy mode in the optical conductivity 
spectrum. 



VII. SUMMARY AND OUTLOOK 

We have performed far-field infrared ellipsometry 
and reflectance measurements on VO2 films across the 
insulator-to-metal transition (IMT) and have determined 
the evolution of the effective optical constants. In ad- 
dition, we have investigated the IMT with scanning 
near-field infrared microscopy and have directly observed 
the percolative phase transition. A combination of the 
far-field and near-field results within the ambit of the 
Bruggeman effective medium theory uncovers enhanced 
effective mass in the incipient metallic islands at the on- 
set of the IMT. This result signifies the pivotal role of 
electronic correlations in the IMT physics. 

We find that the dipole model of the near-field in- 
frared contrast provides a satisfactory explanation of the 
observed amplitude contrast between the metallic pud- 
dles and the insulating regions. The agreement between 
experiment and the dipole model is evident only after 
the near-field amplitude from the metallic regions is nor- 
malized to the near-field amplitude from the insulating 
regions. Future s-SNIM experiments on VO2 with Au 
reference should help better understand the evolution of 
the near-field amplitude as the sample goes through the 
IMT. 

The methods and results reported here set a new 
standard for investigating the conducting state in 
other doped correlated insulators, for example the 
cuprates and manganites, where phase segregation is 
commonplace! 34 | 35|36|37 | 38 | 39 l Ideally, one wants to obtain 
broad-band near-field infrared characterization on the 
nano-scale, without recourse to effective medium theory. 
The technology for this pu rpose is in its infancy and re- 
quires development ! 71 l 72 l 73 l Our work provides an incen- 
tive for such endeavors. 
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IX. APPENDIX A 

The measured ellipsometric coefficients ^f(oj) and A(cj) 
for T = 342 K are compared in Fig. [6^i,b with the results 
from Bruggeman effective medium theory for an inho- 
mogeneous two-phase system. The latter is a fit to the 
experimental data based on the optical constants of the 
monoclinic (Ml) insulator and the rutile metal at T — 
360 K. The filling fractions and depolarization factor used 
for the fits are taken from Table [T] Satisfactory agreement 
is obtained at higher frequencies (10000 - 20000 cm -1 ) 
as can be seen clearly in the plots of the difference be- 
tween the fits and the data displayed in the insets of 
Fig. [6^,b. However, the Bruggeman EMT should pro- 
duce the most accurate results at low frequencies which 
is not the case here. This implies that the optical con- 
stants of the metallic phase in the inhomogeneous IMT 
regime may be different from those of the rutile metal 
at T = 360 K provided the optical constants of the Ml 
insulator do not change in the IMT regime. 



in topography. However, we were unable to study our 
hypothesis that the metallic islands are seeded at grain 
boundaries, defects etc. This is because extensive image 
processing is required in order to identify the relationship 
between the new-born metallic islands and grain bound- 
aries. This analysis can be reliably applied only if the 
topography data are more detailed. Optimized topogra- 
phy images will be required to investigate possible cor- 
relations between the nucleation sites of metallic islands 
and topographic features like grain boundaries. 



XI. APPENDIX C 

The following equations of the extended Drude formal- 
ism ^ describe the procedure for extracting the scat- 
tering rate X/t(lj) and the mass enhancement factor 
m*(u>)/mb from the real and imaginary parts of the com- 
plex optical conductivity of the metallic regions d a {oj)= 
a la (oj) + i(J 2 a{u)- 



X. APPENDIX B 

Near-field amplitude images obtained while cooling the 
V0 2 film through the IMT are shown in Fig. [7] The 
images show that the insulating phase nucleates in the 
metallic host, and then proliferates and finally percolates. 
The IMT occurs about 8 K lower in temperature than on 
heating. Such a hysteresis is expected for a first-order 
phase transition and is seen in FigjlJ). 

In Fig. [8^L,t> we demonstrate that repeated near-field 
imaging at fixed conditions yields unchanged patterns of 
metallic puddles. Imaging takes about 5 minutes (« 0.02 
s per pixel of ~ 15 nm x 15 nm size) and was immediately 
repeated. The reproducible patterns of metallic puddles 
rules out the possibility of slow dynamic fluctuations. If 
dynamic fluctuations occur on shorter time scales, then 
they do not affect the average pattern, especially they do 
not nucleate and stabilize into static metallic puddles. 70 
We hypothesize that the specific patterns that emerge 
are most likely seeded and controlled by strain at the in- 
terface, defects, grain boundaries etc. These nucleation 
centers should lead to a reproducible pattern at a par- 
ticular temperature, a feature that could not however be 
proven in the experimental setup because thermal drift 
due to thermal cycling through the phase transition pre- 
vented the repeated s-SNIM imaging of the same sample 
area. 

We note here that simultaneous topography images 
were obtained along with near-field infrared images dis- 
played in this work. We verified that the near-field in- 
frared contrast between insulating and metallic regions 
in the IMT regime is due to the difference in optical con- 
stants and not due to topographic artifacts or changes 



1 



m* (lu) ^ 



p 

2 



Q"2oM 
M + °2a( W ) 



(7) 
(8) 



The plasma frequency u> p for the metallic states is ob- 
tained from the partial sum rule: 



ai a (u)du> 



(9) 



Here, the upper frequency limit lo c = 13700 cm -1 
is chosen for the rutile metal to ex clude contributions 
from higher- lying optical transitions! 18 ! 19 ! This gives ui p 
— 22000 cm -1 for the rutile metal and this value is nearly 
constant (within ten percent) for the metallic state in the 
nanoscale puddles with fixed uo c = 13700 cm -1 . With 
this nearly constant choice of the plasma frequency for 
all metallic states, the mass enhancement factor for the 
strongly correlated metallic puddles can be directly com- 
pared to the mass enhancement factor of the rutile metal. 

The plasma frequency is related to the carrier density 
(n) and the band mass (m,b) by the following expression: 



47rne 2 



(10) 



We assume that the carrier density does not change 
significantly in the strongly correlated metallic puddles 
compared to the rutile metal. 
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FIG. 1: (color online)(a) Room temperature X-ray diffraction data for the (200)-oriented MI-VO2 film on (1012)-oriented 
sapphire substrate, (b) The resistance of the VO2 film is plotted while heating and cooling across the IMT. (c) Real part of the 
conductivity <ji(ui) of VO2 plotted as a function of frequency for various representative temperatures while heating the sample 
through the IMT. The open circle denotes the isosbestic (equal conductivity) point for all spectra. (Inset) The temperature 
dependence of the real part of the dielectric function ex in the low- frequency limit (u> = 50 cm -1 ), (d) The optical constants 
of VO2 as in (c) but with phonon contributions subtracted. 



FIG. 2: (color online) Images of the near-field scattering amplitude S2 (scale in relative units) over the same 4 /j,m by 4 /im area 
obtained by s-SNIM operating at the infrared frequency to = 930 cm" 1 . The images were obtained for various temperatures 
while heating VO2 through the IMT. The metallic regions depicted in light blue (light gray) and white give higher amplitude 
compared to the insulating phase shown in dark blue color (dark gray). 



FIG. 3: (color online): (a) The distribution of near-field amplitude S2 from images in Fig. [2] for selected temperatures. Note 
the log scale on the vertical axis for clarity of presentation, (b) the distribution of the normalized near-field amplitude for the 
same temperatures as in (a). The median near-field amplitudes for the insulating phase were used for normalizing the near-field 
amplitude in (b) . The distributions plotted in panels (a) and (b) are the number of pixels normalized such that the area under 
each of the histograms is unity. 



FIG. 4: (color online) Optical characterization of the metallic puddles derived from experiment: (a) the real part of the optical 
conductivity ai a (uj) (b) the scattering rate l/r(w) and (c) the optical effective mass normalized by the band value m*(ui) / rrib 
of the metallic regions of VO2 for representative temperatures in the IMT regime. The inset in (c) shows the uo — > limit of 
the mass enhancement factor as a function of temperature (data points between T = 400 K and 550 K are from Ref ll8p . 



FIG. 5: (color online) Dependence of the real part of the optical conductivity ai a {<jj) of the metallic puddles on the depolarization 
factor q, with / fixed at 0.31 for T = 342.6 K. 



FIG. 6: Panels (a) and (b) display the ellipsometric coefficients ^(uj) and A(ui) respectively that were measured for the VO2 
film on a sapphire substrate at T = 342 K (dashed gray curves). The thin black solid lines are the ellipsometric coefficients 
derived from the Bruggeman effective medium theory model as explained in Appendix A. The insets show the difference between 
the fits and the data for the respective ellipsometric coefficients. 



FIG. 7: (color online) Near-field amplitude S2 images (scale in relative units) obtained at u) = 930 cm 1 with decreasing 
temperature through the IMT. 



FIG. 8: (color online) Near-field amplitude S2 images at w = 930 cm 1 at T = 342 K. Image (b) was obtained immediately 
after image (a). 



11 



Fig. 1 




Fig. 4 




Fig. 5 




6 




Fig. 7 




